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A liquid poured from a curved solid surface can separate as a steady jet or sheet, or trickle down 
along the solid surface. It was shown by Duez et al. [Phys. Rev. Lett. 104, 084503 (2010)] that 
surface wettability controls the separation of an inertial flow from a solid surface to an unexpected 
degree, which was further motivated by an inertial-capillary adhesion model. In this paper we extend 
the analysis by a control volume calculation that takes into account the velocity profile within the 
flowing layer, supported by Boundary Integral potential flow simulations, and the detailed capillary 
forces induced by the local curvatures of the sheet. Our analysis captures the appearance of a 
critical Weber number below which no steady separated solutions can be sustained. We investigate 
the dependence of the critical Weber number on the wettability and sharpness of the edge of the 
curved solid, and recover the key experimental trends. 


I. INTRODUCTION 

The so-called ‘teapot effect’ is a daily life phenomenon 
that will be recognized by everyone. When a liquid is 
poured too slowly from the (hydrophilic) nose of a teapot, 
a bottle or a beaker, the liquid has the tendency to run 
down along the underside of the spout, as depicted in 
Fig. [T}r. This remarkable effect already received atten¬ 
tion in the 1950s [l, 2|, and was explained by Keller and 
Vanden-Broeck as a purely hydrodynamic (Bernoulli) 
principle M. Since there is a difference in flow ve¬ 
locity between the top and bottom of the liquid film, the 
pressure is lowest directly above the spout lip; sufficiently 
low that the flow is pushed further down along the convex 
solid. In this situation the liquid film completely ‘trickles’ 
around the solid, and exact potential flow solutions were 
found in which the edge of the solid was treated as per¬ 
fectly sharp ai- After trickling, a possible separation 
from the solid could be induced by gravity. 



FIG. 1. Problem sketch, (a) Water flow trickling around 
the curved (hydrophilic) solid surface, known as the ‘teapot 
effect’; (b) The water repellence of the solid (or the flow speed 
of the liquid) overcomes the trickling behavior. 

However, one can overcome the trickling scenario by in¬ 
creasing the flow speed, such that the liquid jet or sheet 
separates from the solid with an angle, as is sketched 
in Fig. [l]). Simulations of pouring flows have revealed 
complicated dynamics, particularly when the Reynolds 
number of the flow becomes of order unity M- More¬ 


over, Duez et al. [?] found that surface wettability con¬ 
trols the flow separation to an unexpected degree, also at 
large Reynolds numbers. This was also confirmed by re¬ 
cent experiments that aim to develop controlled overflow 
by surface manipulation [§}. Since wettability is deter¬ 
mined on a microscopic scale, this highlights the impor¬ 
tance of multiple length scales in the problem of pouring 
flows. In some sense, the experiments and analysis by 
Duez et al. have revealed a relatively simple framework 
for understanding the various regimes of pouring flow. 
The experiments were performed using an axisymmet- 
ric set-up that is the top part of the ‘liquid bell’ geom¬ 
etry ({Mil}. for which the relevant parameters are the 
flow speed, the wettability of the solid, and the radius of 
curvature at the edge of the impacter. The transition 
between inertia-induced sheet separation and trickling 
appeared at well-defined Weber numbers, whose critical 
values depend on the contact angle and sharpness of the 
solid edge. Using a scaling argument for the horizontal 
momentum balance, Duez et al. predicted a dependency 
on the flow speed of the liquid and the wettability and 
sharpness of the solid surface, capturing the main exper¬ 
imental trends - though surprisingly, the experimental 
trends by Dong et al. [8| suggest a different dependence 
on the curvature of the edge (linear instead of quadratic). 
The modeling approach has so far been limited to a force 
balance in the horizontal (not in the vertical) direction 
and it has remained a challenge to explicitly capture the 
origin of the trickling transition. In addition, it is not 
clear how these results relate to the analytical potential 
flow solutions ii , which excluded the effect of surface 
tension and finite edge-curvature. 

In this paper, we will perform a control volume anal¬ 
ysis for inertial pouring flows over a solid edge of finite 
curvature, by taking into account the capillary forces in¬ 
duced by the shape of the meniscus and the velocity pro¬ 
file within the liquid. By releasing some geometric con¬ 
straints imposed in Q, this will allow for a force balance 
in the horizontal and vertical direction. The model is 
solved for varying flow parameters and geometric param- 
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eters, and we show how this indeed leads to the appear¬ 
ance of a critical Weber number for the trickling transi¬ 
tion. In Sec. m we will introduce the relevant parameters 
in the problem and the basic assumptions. This results 
into a set of coupled equations that provide a predic¬ 
tion for the separation angle. Our main findings will 
be presented in Sec. uni where we give both numerical 
and asymptotic predictions for the critical Weber num¬ 
ber. The results will be summarized and compared to 
experiments in Sec. IIVI 


II. MODEL 

A. Definitions and assumptions 

A sketch of the problem and the relevant parameters 
is given in Fig. HJi. Like Ref. [?], we treat the flow as 
two-dimensional and define a horizontal x direction, and 
a vertical y direction. We assume a steady, laminar, ir- 
rotational flow with high Reynolds number. In that case 
the flow can be considered uniform both at the inflow 
above the solid and within the separated sheet, with ve¬ 
locity U and the film thickness is h. The edge of the 
solid has a circular shape, characterized by the radius 
of curvature 77 . We assume the sheet separates at an 
angle a with respect to the horizontal direction, and an 
important part of the analysis is to determine this an¬ 
gle. The sheet separates from the solid at the position 
angle /? defined with respect to the vertical axis. Note 
that, in general, a and we treat these angles as in¬ 
dependent parameters. At the separation point, a small 
capillary meniscus is formed, which has a radius of cur¬ 
vature r m , which is set by the Yo ung -Laplace pressure 
difference over the free surface [l2|, El- Locally, the cir¬ 
cle formed by the meniscus crosses the 77 -circle with the 
contact angle 6q, which is how the surface wettability en¬ 
ters the analysis. Note that r m is typically much smaller 
than 77 (Fig. [5] not drawn to scale). 

Assuming inviscid flow, the introduced parameters can 
be expressed by three dimensionless control parameters, 
which are the Weber number We = pU 2 h/ 7 , where p and 
7 are the density and liquid-gas surface tension of the 
liquid, respectively, the ratio 77 = 77 // 1 , and the contact 
angle Oq. For characteristic values around U ~ 1 m/s 
and h ~ 0.1 — 1 mm, the Weber number for water is of 
order 1 to 10. The corresponding Reynolds numbers are 
about 10 2 - 10 3 , so we can indeed assume laminar flow 
and neglect viscosity. 

Clearly, the formulation presented here closely follows 
Duez et al. Q, though some notable differences appear. 
Importantly, we release the geometric constraint that the 
tangent of the sheet (with angle a with respect to the 
vertical direction), is also a tangent line to the edge of 
the solid (i.e. the circle of radius 77 ). In this manner 
we can treat a and fi as two independent parameters, 
and the momentum balance can be maintained in both x 
and y directions. Other minor differences appear below, 


when estimating the forces acting on the control volume. 

We will solve the deflection angle a of the liquid sheet 
as the result of the hydrodynamic and capillary forces 
using linear momentum conservation. The mechanism 
responsible for the liquid sheet bending is that the veloc¬ 
ity profile over the curved solid is not uniform, contrarily 
to the inflow and the outflow of the control volume indi¬ 
cated in Fig. [2] This velocity profile implies a pressure 
difference across the film, which induces a hydrodynamic 
force, Fd , exerted by the solid on the liquid, pulling the 
liquid along the curved edg e. This is sometimes referred 
to as the Coanda effect [14| . The resultant capillary force 
in the small meniscus is denoted as F men , while the effect 
of surface tension over the top of the curved sheet yields 
F top . Importantly, we can define the CV in two ways, 
indicated by Fig. [2Ji and b, respectively. In Fig. [5K the 
Control Surface (CS), is located just inside the liquid; in 
Fig. [2 ]d the CS is located just outside the liquid. In the 
first case, the capillary force can be evaluated by integrat¬ 
ing the Young-Laplace pressure over the control surfaces; 
in the second case, there is no normal force working on 
these surfaces, but the surface tension acts parallel on 
the edge of the liquid domain. Both points of view of 
course result into the same final equations and the resul¬ 
tant forces are indicated by the arrows in Fig. [2 Jd. The 
resulting momentum balance for the x direction reads 


pU h (cos ex 1) — Fd tX + F rnen ^ x + Ft 0 p' X , (1) 
while for the y direction we have: 


pU h sin a — ^d,y F F men y T Ft op y. (2) 

Here we used that the Control Volume (CV) contains 
the full bending of the stream, such that the left and 
right boundary of the CV are located where the flow is 
uniform. 

It is instructive to consider the momentum balance in 
the context of the exact solutions by Keller & Vanden- 
Broeck M- These are obtained by treating the solid 
as a perfectly sharp edge in the absence of capillary ef¬ 
fects. This effectively corresponds to We = 00 , 77 = 0, 
while the angles 6q and fi are not defined in this limit. It 
was found that analytical solutions exist for each value 
of a; one could use the momentum balance (HU with¬ 
out the capillary forces ( F men = F top = 0) to determine 
the dynamical force Fd, for each value of a. Importantly, 
however, this analysis does not lead to a selection of a. 
A selection of a does appear at finite values of We and 
77 , one thus needs explicit expressions for all terms in 
Sys. (11I2D . as will be developed below. In Sec. IIIB1 we 
investigate the velocity profile and the pressure distribu¬ 
tion in the liquid, from which we compute the various 
forces in Sec. Ill Cl The resulting set of equations will be 
presented in Sec. Ill Dl 
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FIG. 2. Flow around the edge of a solid and two different definitions of the Control Volume (CV), indicated by the gray area, 
(a) Definition sketch and the Control Surface (CS) located just inside the liquid. A liquid film with thickness h, which has a 
uniform flow velocity U, bends around a solid surface with radius of curvature ri. do is the solid-liquid contact angle, a is the 
global deflection angle of the separated liquid sheet with respect to the horizontal (x) direction, and f3 is the angular width of 
the wetted fraction of the solid, a and /? are a priori independent. r m is the radius of curvature of the circle formed by the 
meniscus. The figure is not drawn to scale, typically r m -C ry. The value of the marked angle * is n — 8o—/3+a. (b) The force 
balance and the CS located just outside the liquid. The forces, indicated by the arrows, respectively are the capillary forces 
acting in the interface (denoted by 7 ), and the resultant force F,i induced by the pressure difference over the sheet ( Pa — Pb )■ 


B. Flow profiles and pressures 

To compute the pressure distribution in the liquid we 
require the velocity profile inside the flowing liquid, for 
a given geometry characterized by 77 . Here, we focus on 
the case q > 1 for which the flow will evolve towards 
concentric ‘circular’ streamlines around the circular edge 
of the solid 0 - Under the assumption of potential flow, 
this corresponds to a free vortex with a tangential flow 
velocity u ~ 1/r. Such a profile is quite different from 
the corner solutions for f, <C 1 , and the resulting pressure 
distribution is expected to be quite different. 

We therefore verified the 1/r profiles for an experi¬ 
mentally relevant case, f,; = 4, We = 55, using potential 
flow simulations using an axisymmetric Boundary Inte¬ 
gral (BI) routine (lT Il8l | . We solve the Laplace equation 
V 2 ip = 0 for the flow potential tp in the domain indi¬ 
cated in Fig. [3]a, containing a small inlet region before 
the circular bend, ending with a separated sheet. In the 
simulations the contact line is pinned at a fixed position 
angle /?. Note that the BI simulations are only used for 
the confirmation of the presumed velocity profile - we 
have not succeeded in creating perfectly steady sheets in 
the simulations, except for trivial solutions. 

The result is shown in Fig. [3 ]d, where we plot the nor¬ 
malized velocity profile across the film for various loca¬ 
tions 7 / along the curved edge (defined in Fig. [3Jr). The 
initially uniform inflow indeed rapidly evolves towards a 
1/r profile, indicated as the dash-dotted line. The en¬ 
trance/exit effects are quite small: for the example of 
Fig. [3] it is not more than a few degrees at the beginning 
of the bend, and about 10 degrees at the bottom part 
where the liquid separates from the solid. This shows 
that a 1/r profile for f \ 1 is a valid approximation. 

Thus, we base our further calculations on the velocity 


profile u = A/r , were the constant A is determined from 
mass conservation, 


Uh = 



( 3 ) 


so that 


u = U 



( 4 ) 


The flow velocity directly above the solid surface is larger 
than at the top of the liquid sheet, so the pressure dif¬ 
ference P A — Pb (cf. Fig. (2Ji) is indeed positive from the 
steady Bernoulli equation: 


A p = p a -p b = -p (u\;L r . - u\ 2 r=r . +h ) 


= P U 2 


(1 + h) 2 (ln(! + i)) 


( 5 ) 


Computing P A (gauge) using the Young-Laplace pressure 
difference as 


Pa 


1 

h 


1 

1 + fi 


( 6 ) 


we obtain 


Pb 


7 

h 


1 

1 + f% 


pU 2 


QinY 

fi 


( 7 ) 
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FIG. 3. Potential flow simulation of the velocity profile in the flowing layer, (a) Sketch of the simulation domain. The red 
lines are undeformable solid boundaries, while blue edges are deformable and exhibit capillary pressure. The dots are nodes of 
the boundary integral simulations. Parameters are taken f, = 4, We = 55. The simulated setting is in fact axisymmetric, like 
the experiment using an impacter of radius Rim P 0 ; here we adapted a typical experimental value Rim P /h = 26. (b) Velocity 
profiles (normalized by the average U a ve), measured across the liquid film for different locations around the curved solid. The 
velocity profile quickly evolves towards the expected 1 jr profile (dash-dotted line). 


Here we introduced a dimensionless geometrical factor 


where the integration variable ip is the angle with respect 
to the vertical (see Fig. [3ji). Similarly, we find for F ( i y \ 


G(fi) 


_ \ + h _ 

nO + hf (in(l + i )) 27 


( 8 ) 


Fd.y ~ 


PsdAu 



PsTi COS ip dip 


which encodes all relevant information of the velocity pro¬ 
file. 

Using that Pb = —7 /r m , we also extract an expression 
for r m : 


^ m 

h 


_ 1 + fj 

We (1 + J-) g(n) 


- 1 


(9) 


pu 2 hg(fi) -7 



sin/?. 


( 11 ) 


The capillary forces are induced by the Young-Laplace 
pressures over the free surface of the meniscus ( F m en,x 
and F mentV ), and over the top free surface of the film 
(F toPtX and F topy ). For the capillary force induced by 
the curvature of the top of the sheet, projected in the x 
direction, we find 


Thus, the size of the meniscus r m adapts to accommodate 
the low pressure Pb ■ Note that We is typically large, and 
hence r m <C 77 . 


C. Calculating the forces 


Ftop,x — 


nOt 

/ Pa (r% + h) si 
Jo 

fG 


sin ip dip 


Ji + h 
= 7 [cos a — 1 ], 


(77 + h) sin ip dip 


( 12 ) 


We can now proceed to the evaluation of the several 
terms in the momentum balances m by integrating the 
local pressures along the different sections of the CS of 
Fig. [5Ji, which is located just inside the liquid. Along 
the solid, we find for the force induced by the pressure 
difference projected in x direction, F,i.x : 


Fd,x = 


P B dA we t = / PbU sin ipdip 

- x J 0 


pU 2 hg(ri) - 7 


1 + fi 


(cos P — 1), (10) 


and 


Fto P ,y = -7 sin a. (13) 

Note that the expressions for F toPtX and F top! y are more 
easily interpreted by considering the explained equivalent 
picture of the force balance in Fig.(2J:>: the separate terms 
are the x and y projections of the capillary forces along 
the free surfaces of the CV. A similar observation holds 
for the capillary forces in the meniscus, which can be ob¬ 
tained by integrating the pressure Pb over the meniscus 
circle with radius of curvature r m 





















5 


-(7T-/3-0O) 


P B rm sin ip dill 


7 


r m sin ipdi/j 


= —7 [cos (—(n — /3 — 0o)) — cos a] 
= 7 [cos(/? + 9 0 ) + cos a], 


(14) 


and 


Fmen,y = 7 [- sin(/3 + 6 0 ) - sin a]. (15) 


Here we used the fact that the meniscus angle, * in Fig.[2j 
is equal to 7 r — 9$ — /3 + a, as follows from the geometry. 


D. Resulting system of equations 


The momentum balance m combined with the com¬ 
puted forces finally yield the key equations of the model: 


We — Wet/(fj) 


1 

(1 + n). 


[2 — We] cos a + 


[2 — We] sin a + 


with We = pU 2 h/"f, as previously defined. This system 
should be seen as equations for the angles a and /3, which 
can be solved for given values of We, iy, and do- Note 
that by imposing a contact angle 9q, we by definition 
consider only separated sheets. 

Before we proceed to the numerical solutions of (11611711 
in Sec. IIIII we can already analyze a few interesting lim¬ 
its. The limit We —> oo gives as only solution a = /3 = 0, 
corresponding to the perfect horizontal sheet, for every 
value of do- Interestingly, the same holds for do = 180°, 
for which there is no capillary adhesion: if the surface 
is superhydrophobic, the jet/sheet will be perfectly hor¬ 
izontal for any We. A third interesting case is fj —>• oo, 
for which G(fi) —>• 1 and no separated solution exists, 
unless do = 180°. These trends all agree with the exper¬ 
iments in Ref. 0. A final special case is We = 2. In 
that case, the a dependence completely drops out of the 
equations, and we are left with two equations for a single 
unknown, /3. This has no solution unless do = 180°. We 
can interpret We = 2 as a minimum flow speed needed 
for a non-retracting sheet, as this indeed coincides with 
the Taylor-Culick velocity [l3, 20]. 


III. RESULTS 

A. Solutions 

We now analyze the solutions of the momentum bal¬ 
ance (1161171) . In Fig. 0Jr we report the separation angle 
a as a function of We for fj = 4. The various curves 
correspond to different wettabilities, with d 0 increasing 
along the arrow. It is found that solutions only exist 
above a critical value of the Weber number, We c , which 


We<7(fj) 

We^(fj) 


Tj 

(i + u). 


cos/? + cos(/? + do) = 0 ; 


Tj 

(1 + n)_ 


sin/3 + sin(/3 + d 0 ) = 0, 


(16) 

(17) 


we identify as the threshold for the trickling transition. 
The critical point is found to coincide with a = 180°. 
Above We c , the momentum balance admits two possible 
solutions. However, solutions for a larger than 180° are 
not physical in the sense that the liquid would cross the 
solid and we focus on the lower solution branch. As ex¬ 
pected, the deflection angle a increases when the fluid’s 
inertia is reduced, i.e. as the Weber number is decreased. 
The sketches in Fig. 0)D-d further illustrate this effect. 
It should be emphasized that the deflection angle varies 
rapidly with We for a beyond 90°, i.e. when the crit¬ 
ical point is approached. The inset of Fig. |4ji shows a 
zoom around the critical point for both angles, a (blue 
solid line) and /3 (red dashed line) for 9q = 90°. The 
two angles always take similar values, with a maximum 
difference of about 20°. The global minimum of /? is also 
reached at We c , but has a value slightly below 180° (equal 
to 180° if 9 0 = 0 °). 

The critical Weber number is thus found to depend on 
the wettability of the solid. Figure [5] presents We c as a 
function 1+cos 9q for several values of f, (solid lines). The 
dependence is approximately linear, with a slope that 
increases for increasing fj. The linearity no longer applies 
when 1 + cos 9o is small, i.e. the superhydrophobic limit, 
for which all curves converge to We c = 2. In this case 
there is no adhesion by the solid and the minimal velocity 
is given by the Taylor-Culick velocity of a free sheet. 

Figure [ 6 ] shows We c as a function of the radius of cur¬ 
vature of the solid fj, for several values of 9q (solid lines). 
For large fj we again find a linear trend, reflecting that 
it is more difficult to separate the sheet when the edge 
of the solid is not sharp. Note that the model is only 
valid for fj > 1 , owing to the assumptions of the velocity 
field that leads to (| 8 |). It is still interesting to report the 
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FIG. 4. (a) q vs We = pU 2 h /7 for ft — 4 and different 

contact angles. Following the arrow, do increase from 0° to 
150° in steps of 30°. Solutions only exist for We larger than a 
critical Weber number We c . Inset: a detailed view of a and /3 
as a function of We around We c for 8q — 90°. (b-d) Resulting 
flow contours, corresponding to the marked dots in panel (a), 
showing the dependence of the separation angle on the Weber 
number for do = 150°. 


model prediction for r, < 1: there is even a divergence 
that appears when Q takes the value 1 (Fig. [6l dotted 
lines). In the Discussion section we briefly comment on 
the limit of small fj. 

B. Asymptotic expansion 

It is instructive to attempt an approximate solution 
for the critical Weber number, based on the observation 
that at the critical point a = 180° and /3 = 180° — e, 
where e turns out to be small (< 20°). Hence, we ex¬ 
pand Sys. (1161171) around up to first order in e, replacing 
sin e ~ e and cos e ~ 1. Using sin (0 O — e) = sin d 0 cos e — 
cos 9 0 sin e and cos (do — e) = cos do cos e — sin d 0 sin e, we 
find for e (in radians): 


G(h )We + cos do - j^rr' 

and 

We)? + AWe c + B « 0, (19) 

with 



FIG. 5. Critical Weber number We c as a function of 1 + cos do 
for n = 0.5, 1, and 2 mm, and h = 0.25 mm. Solid lines are 
numerical solutions of the momentum balance; dashed lines 
correspond to the approximation 09 - 



fig. 6 . Critical Weber number We c as a function of r, = n/h, 
h = 0.25 mm, for do = 10°, 90°, 120°, and 175°. Solid lines are 
numerical solutions of the momentum balance; dashed lines 
correspond to the approximation ®. The model is not valid 
below ri ~ 1 , where the predictions are given by a dotted line. 


J 


2 (t^-) + cos do (3 S(fi) - 2) + 2Q(r i ) (l - + G{h) (^ff) 

A ~ -2G(fi) + 2 G(n) 2 


( 20 ) 


and 


B = 



2 cos d 0 


fi(fj-i) 

(i +U? 


(t+w) cos °o + 1 


-2G(r t ) + 2G(n) 2 


( 21 ) 
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This gives a quadratic equation for We c , which is su¬ 
perimposed as dashed lines in Figs. [5] and [G] Indeed, this 
approximate solution gives a very good description of the 
full solutions. The difference is largest for 9q ~ 90°, where 
the largest values for e are encountered. 

We could even further simplify (H9l) when f1. G{fi) 
can then be approximated by 1—1/ (2fj)), and this yields 

We c ~ fj (1 + cos 9q) . (22) 

This predicts linear behavior of the critical Weber num¬ 
ber with respect to both the aspect ratio fj and to 
1 + cos 0o, which is consistent with our results in Figs. [5] 
and [6] 


IV. DISCUSSION 

We have performed a control volume analysis of in¬ 
ertial pouring flows, taking into account the fluid iner¬ 
tia, surface tension, and the curvature of the edge of the 
solid. The analysis explicitly recovers that steady states 
can only exist above a critical Weber number We c , and 
thus captures the transition to trickling. The work also 
recovers the experimental trends that We c decreases as 
the solid is sharper and more hydrophobic, and identify 
a lower bound We c > 2. Here we wish to conclude by 
making a direct comparison to experiments [?, 8j. 

Figure [7}i reports the separation angle a against 
We/We c for three different contact angles, 9 0 = 10°, 
115°, and 175°. The value of fj = 4 was adapted from 
the experiment, and the symbols represent data from 
Ref. Q. It is clear that the model captures the exper¬ 
imental trends, but is limited in terms of quantitative 
prediction. In particular, the change of angle a close to 
the critical point is underpredicted by the model. Also, 
the experiments at large Weber number exhibit a sheet 
deflection of the order of 10°, which points to either dis¬ 
sipation in the fluid or an influence of the axisymmetric 
setup - both of which are not taken into account in the 
model. 

In Fig. [ft, we show the critical Weber number We c as 
a function of 1 + cosd, for three different fj. The model 
nicely captures the linear dependence We c ~ (1 + cos 9 o) 
that is observed in experiment. In addition, the data in 


Fig. [Tji are consistent with the predicted lower bound on 
the critical Weber number of We c = 2, set by the Taylor- 
Culick velocity. The model underestimates the depen¬ 
dence on fi : we predict a linear increase of We c with f j, 
while the experimental data are better described by r~, 
as was also proposed by a scaling argument 0|. However, 
the recent data in Ref. § are best described by a linear 
dependence on fj, even though their experimental setup 
is in principle comparable. At present, there is thus some 
uncertainty on how the sharpness of the edge affects the 
trickling transition. We do note that the model predic¬ 
tions shown here are without adjustable parameters and 
the discrepancy with, and between, various experiments 
is always less than a factor of two. 

The momentum balance presented in this paper has 
the merit that it provides a general framework for future 
investigations of pouring flows. Namely, the capillary 
forces are described correctly here - apart perhaps from 
the question of whether 9q can be interpreted as the equi¬ 
librium contact angle [2ll . [22f - and the main assumption 
is on the estimate of the hydrodynamic force Fd . As 
mentioned already, we only estimated this force in the 
case of large fj, but IMP can in fact be used to estimate 
Fd based on experimental data. The black triangles in 
Fig. [7]) correspond to experiments with fj = 0.12, for 
which We c exhibits almost no dependence on wettabil¬ 
ity and is close to the lower bound We c = 2. The mo¬ 
mentum balance suggests that F ( j —> 0 in this limit, im¬ 
plying a vanishing hydrodynamic retention around sharp 
edges. Interestingly, this is a strong departure from the 
exact potential flow solutions around a perfectly sharp 
edge [1, H, which do not include capillary adhesion, for 
which trickling completely relies on a nonzero Fd■ Fu¬ 
ture work should further reveal how the presence of a 
meniscus on a sharp edge influences trickling in inertial 
pouring flows. 
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